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ABSTRACT 

Solar p-mode frequencies vary with solar activity. It is important to take this into account when comparing the frequencies observed 
from epochs that span different regions of the solar cycle. We present details of how to correct observed p-mode frequencies for 
the effects of the solar cycle. We describe three types of correction. The first allows mode frequencies to be corrected to a nominal 
activity level, such as the canonical quiet-Sun level. The second accounts for the effect on the observed mode frequencies, powers, 
and damping rates of the continually varying solar cycle and is pertinent to frequencies obtained from very long data sets. The 
third corrects for Sun-as-a-star observations not seeing all components of the modes. Suitable combinations of the three correction 
procedures allow the frequencies obtained from different sets of data to be compared and enable activity-independent inversions of 
the solar interior. As an example of how to apply the corrections we describe those used to produce a set of definitive Sun-as-a-star 
frequencies. 
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1. Introduction 

The Birmingham Solar Oscillations Network (BiSON) observes 
Sun-as-a-star Doppler velocities. BiSON data can, therefore, be 
used to determine the frequencies of low-degree (low-/) oscil- 
lations. Low-Z oscillations are particularly important as they are 
the only p-modes that penetrate the solar core. 

Accurate and precise solar p-mode frequencies are often de- 
termined from time series that contain many years of data. One 
drawback to examining such long data sets is that, because of 
the 11 -yr solar cycle, they must span a wide range of different 
activity levels. Since solar p-mode frequencies, lifetimes, pow- 
ers, an d peak asym metries all vary systematically with solar ac- 
tivity dHowd 12008, and references therein) care must be taken 
when analysing the mode parameters obtained from long sets 
of data. For example, between solar maximum and minimum 
the frequency of a low-/ mode at approximately 3000 yuHz can 
change by as much as 1 //Hz, an amount that is far greater than 
the errors associated with fr equency estimates. 

iBroomhall et al.l (120091) has recently published a list of 
definitive low-Z frequencies that were obtained from 8640 d 
of BiSON data. The authors list the raw frequencies, which 
were determined directly fr om the data by fi tting profiles to the 
frequency-power spectrum. IBroomhall et al. also quote frequen- 
cies that have been corrected for various solar cycle effects. This 
paper describes, in detail, how these corrections were made in 
a manner that will allow the reader to reproduce the corrections 
with their own data. This will allow comparisons between the 
frequencies obtained from data sets that were observed at differ- 
ent epochs and consequently different activity levels. 

We begin, in Sect. [2] by describing the well-known linear so- 
lar cycle correction that can be used to correct mode frequencies 
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to a nominal activity level. Then, in Sect. [3] we describe how to 
determine the 'devil-in-the-detail' correction that accounts for 
cross-talk between the variations of different mode parameters 
with the solar cycle, and the distribution of activity levels over 
the period of observations. Finally, Sect.|4]describes a correction 
that is necessary when using Sun-as-a-star observations because 
not all mode components are detectable. 



2. Linear solar cycle correction 

The raw frequencies determined for two sets of data will be dif- 
ferent if they were observed at different epochs. To allow fre- 
quencies from different periods of time to be compared, a so- 
lar cycle correction must be performed. This correction also 
allows the frequencies that would have been observed at the 
canonical quiet Sun level to be determined. The correction is 
based on the assumption that variations in global activity in- 
dices can be used as proxies for low-Z frequency shifts. We 
also assume that the correction can be parameterised as a lin- 
ear functi on of the chosen activity measure. For the 10.7-cm ra- 
dio flux (Tapping & Detracev 1990), these assumptions are ro- 
bust ( Chaplin et al. 2004c) at the level of precision of the data. 
IChaplin et alJ (120071) find that the NOAA Mg II H core-to-wing 
ratio dViereck et aL |2001|) is predominantly a better proxy for 
making solar cycle corrections. However, we found that the fill 
of the Mg II data was too sparse to use here. 

Let v n j be the set of fitted frequencies determined using data 
that were observed during a time when the mean activity level 
was (A(t)). For example, the mean 10.7 -cm radio flux observed 
in the 8640-d time series examined by Broom hall et al.1 {2009) 
was (A(f)> = 118x 10 _22 Wnr 2 Hz _1 . (For the remainder of this 
paper we refer to 10~ 22 Wm~ 2 Hz _1 as radio flux units, RFU.) 
This value of (A(f)) has been calculated taking into account 
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Table 1 . Values of gi used to make the linear solar cycle correc- 
tions in IBroomhall et all d2009t) . 

I gi (xlQ-yHzRFU-') 

2.3 ±0.1 

1 3.1 ±0.1 

2 3.1 ±0.1 
2 3.3 ±0.2 
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Fig. 1. Linear corrections that were applied to the frequencies 
observed in the 8640-d time series, where the average activity 
level was 118 RFU, to obtain frequencies correct for the canoni- 
cal quiet-Sun level of 64 RFU. 



any gaps that were present in the BiSON time series (e.g. due 
to inclement weather and very occasionally instrumental prob- 
lems). Let ( A(t)) c be the acti vity level we wish to correct to. 
For example, iBroomhall et alj quote frequencies that have been 
corrected to the canonical quiet-Sun value of the radio flux, 
which is fixed, from historic al observations, at (A(t)) c = 64 RFU 
(Tapping & Detracey] |1990l) . The linear solar cycle correction, 
6v„j, is given by 



Sv nJ = -g,T(v) [<A(0> - <A(f)>J . 



(1) 



where gi are /-dependent factors that calibrate the size of the 
shift and were determined in the manner described below. The 
T(y) in equationQJis a function that allows for the dependence of 
the shift on mode frequency. We have used the same freq uency 
dependence as was determined by IChaplin et al.l d2004blcl) . The 
calculated 6v„j can be added to the set of frequencies, v„j, to 
obtain frequencies that are corrected to a mean activity level of 
<A(f)> c . 

To determine gi the 8640-d set of BiSON data was split into 
20 contiguous 432 d segments. To uncover the dependence of the 
solar-cycle frequency shifts on the 10.7-cm radio flux mode fre- 
quencies were deter mined for each of these 20 time series in the 
manner described by Chaplin ^t alj (12004c). The gradient of the 
linear relationship between activity and frequency varies signif- 
icantly with I because of the spatial dependen ce of the surface 
activit y. Table [TJ gives the values of gi used in IBroomhall et al.l 
(2009). Figure [TJ shows the linear frequency shifts, <5v„j, that 
were applied to the 8640-d BiSON data to correct the frequen- 
cies to a mean activity level of 64 RFU. 

The errors on 6v n j are dominated by the uncertainties asso- 
ciated with g/. These errors need to be propagated through to the 
corrected frequencies and so the corrected uncertainties are of 
the order of 10 per cent larger than those associated with the raw 
fitted frequencies, v„j. 

It is possible that the true relationship between the observed 
frequency shifts and solar activity is quadratic or even of higher 
order. However, we do not find a higher-order relationship when 
we analyse the BiSON frequency shifts i.e., a linear relationship 



provides the best fit, given the precision in the BiSON frequency 
shifts. In order to ensure our correction procedure is internally 
consistent, we use BiSON shifts and a linear calibration. 

An alternative approach is to use se parate correct i ons for 
the rising and falling phases of the cycle. IChaplin et all j2004c) 
show that there is a significant difference in the slope of the lin- 
ear relationship between the rising and fa lling phases of the so- 
lar cycle. More recently. IJain et al] d2009l) found, when studying 
modes of intermediate degree, that different slopes are required 
to accurately describe the observed changes in frequency when 
the cycle is split into different phases. Such an approach would 
be worthwhile if the data for which the solar cycle corrections 
are being made span only a rising or a falling p hase. 

Ho wever, the time series considered by IBroomhall et al] 
(2009) spans 8640 d, including two rising phases and two falling 
phases. To test whether it was advantageous to fit the phases sep- 
arately we have split the 432 d time series into two sets, one for 
the rising phases and one for the falling phases. A linear fit be- 
tween frequency shift and activity was performed for each set. A 
weighted mean of the linear corrections for the two phases was 
then calculated, where the weights corresponded to the propor- 
tion of the 8640 d spent in each phase. We found that for the- 
oretical modes at 3000 //Hz with / in the range 0-3 the two- 
phase frequency correction was in good agreement (within lcr) 
of the corrections found when a single linear fit was applied. 
Furthermore, the error bars associated with the two different cor- 
rection approaches were similar in size. In other words, when 
considering a long time series that spans both rising and falling 
phases of the cycle it is more than adequate to use a single linear 
fit because the differences between the rising and falling phases 
are averaged. 

3. Devil-in-the-detail and activity distribution 
correction 

3.1. A description of the devil-in-the-detail effect 

The linear solar cycle correction assumes that the observed mode 
frequencies correspond to an unweighted av e rage o f the time- 
varying frequencies. However, Chapli net al] (I2008I) found that 
when analysing long data sets this is not the case. Long time 
series span, at the very least, a sizeable fraction of the 1 1-yr so- 
lar activity cycle. The data examined by Broom hall et al] (J2009) 
is a good example of such a time series as it sp ans more than 
two complete solar cycles. IChaplin et al] d2008l) find that the 
frequencies observed in long time series are biased by a cross- 
talk effect that has its origins in the simultaneous variations of 
mod e frequencies , powers, and linewidths over the solar cy- 
cle. IChaplin et al] call this bias to the observed frequencies the 
'devil-in-the-detail' effect and we now describe in detail why 
this effect is observed. 

Consider a simplistic scenario where only the mode frequen- 
cies change with time. The observed frequency of a mode is then 
determined by the activity levels observed over the total period 
of observations and the unshifted frequency of the mode in ques- 
tion. Now consider the fictitious scenario where we have lOOd 
of data and on 99 of the days the solar activity level is constant. 
Then on the 100 th day the activity level doubles. 

The final observed profile of a mode is well represented by 
the sum of the mode profiles observed on each of the lOOd, nor- 
malised by the length of the observations. For the first 99 d of 
our fictitious time series the same profile is observed. However, 
on the last day the peak frequency of the mode profile will be 
shifted by an amount, given by equation [TJ which is dependent 
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Fig. 2. Distribution of t he 10.7-cm radio flux o ver the 8640-d 
time series examined by Broomhall et al. (2009). The solid line 
shows the flux distribution and the dashed line indicates the 
mean flux of the complete time series once the BiSON window 
function has been accounted for (i.e. 118 RFU). 



on the mode's degree and frequency. The influence of the final 
day on the peak frequency of the integrated mode profile will 
depend on the size of the frequency shift experienced by the 
mode and the mode's linewidth. If the frequency shift is signifi- 
cantly smaller than the mode linewidth the final-day profile will 
be approximately equal to the unweighted average of the daily 
mode frequencies. However, if the frequency shift of the mode 
is large compared to its linewidth the peak of the profile for the 
final day will be situated in the wing of the integrated mode pro- 
file. Consequently, the integrated profile's peak frequency will 
not equal the unweighted average of the daily mode frequencies. 
Effectively, the final peak frequency of the integrated profile is 
a weighted average of the time-varying frequencies, where the 
weights are a function of the frequency shift observed on each 
day and the width of the mode profile. 

The above example is extreme and the described effect is 
not important if the activity levels observed during the period of 
observations are uniformly distributed. However, as can be seen 
from Fig.[2]t he activity distribution o ver the 8640-d observations 
examined bv lBroom hall et al.] d2009l) is heavily skewed towards 
low activity levels. At low activity levels the observed frequen- 
cies are at a minimum. Therefore the frequency of the integrated 
profile will be smaller than that expected from an unweighted 
average of the instantaneous frequencies. 

This is not the only effect that influences the final observed 
frequency. As mentioned, it is simplistic to assume that only the 
mode frequencies change with activity. In fact the heights and 
widths of the mode profiles also vary throughout the solar cycle 
and this means that the data observed at different epochs have 
varying contributions to the final peak profile. For typical low- 
/ modes, as the surface activity increases the mode frequencies 
and widths increase but the mode heights decrease. Data samples 
from times when the mode is more prominent (i.e. at low activity 
levels) carry proportionately larger weights in the final peak pro- 
file than times when the mode is less prominent (i.e. at times of 
high activity). Therefore, data samples from times when the fre- 
quencies are lower carry proportionately larger weights in deter- 
mining the final peak profile and so the final observed frequency 
is biased towards a lower frequency than would be observed in 
an unweighted mean. The size of the effect is dependent on the 
comparative time variations in frequency and power. 

We now describe how the magnitudes of both of these ef- 
fects were determined for the 8640-d time series examined by 
iBroomhall et"ai1 d2009l) 



Fig. 3. Corrections applied to the 8640-d frequencies that oc- 
cur due to the skewed distribution of the observed activity level 
i.e. the devil-in-the-detail corrections calculated when the mode 
heights and widths are kept constant with time. 

3.2. How the devil-in-the-detail correction was determined 

Chap lin et alJ (120081) find that the final mode profile observed in 
a long data set is well described by the integral over time of the 
instantaneous Lorentzian profiles sampled at time, f, which have 
the appropriate frequency, power, and damping rate for that time 
and activity level, and so 



{P(V)} = i c T j^l 

T J f=0 !+£(/) 



dr. 



where 



2[v-v(Q] 
A(0 ' 



(2) 



(3) 



Here, H(t) is the instantaneous peak height of the profile, T is the 
total observing time, v(f) is the frequency of the mode at time, f, 
and A(f) is the instantaneous width of the mode profile. 

To determine the size of the frequency bias that occurs be- 
cause of the distribution of activity levels we consider the sim- 
plistic case of a single mode whose peak height, H(t), and width, 
A(f), remain constant with time. However, we take the peak fre- 
quency of the profile, v(/), to be a linear function of the activity 
at time, t . 

For each day in the 8640-d time series we have produced 
a mode profile whose frequency was determined by the 10.7- 
cm radio flux observed on that day i.e. we use equation Q] to 
determine the frequency of the mode at time, t . We then summed 
the individual day Lorentzians to produce a single profile that 
was representative of the profile observed over a long time series. 
We found that the frequency of the final profile was not equal to 
the unweighted average of its time-varying frequency. Figure [3] 
shows the frequency correction applied to correct the observed 
frequency to the unweighted average frequency when H(t) and 
A(f) were constant. As can be seen the correction is significant 
and structured. 

To understand this effect in more detail we consider the fol- 
lowing ratio 



R = 



An ,/ ' 



(4) 



where dv n j was determined using equationQ]for (A{t))-(A{t)) c = 
(1 18-64) RFU = 54 RFU. Figure|]shows R plotted as a function 
of the mode frequency. Comparison with Fig.|3]shows that R has 
a similar structure to the observed frequency shift. Hence, the 
larger the shift experienced relative to the width of a mode the 
greater is the difference in the integrated profile frequency and 
the unweighted-mean frequency. 
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Fig. 4. Variation in R with mode frequency. 
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Fig. 5. Left panel: Model of how the percentage change in the 
mode profile width, A n j, changes with frequency. Right panel: 
Model of how the percentage change in the mode profile height, 
H, changes with frequency. Both panels represent the maximum 
change expected and so are applicable when the activity level is 
at a maximum. 



In reality the above example is simplistic as the heights and 
widths of the mode profiles also vary with the solar cycle. To 
determine the size of the frequency shift due to both the distri- 
bution of the activity and the cross-talk between the cycle vary- 
ing mode parameters we systematically altered the mode profile 
widths and heights with activity. 

Between solar minimum and maximum mode widths are ob- 
served to increase by up to 17.5 per cent a nd mode heights are 
observed to decreas e by twice this amount (jChaplin et alfe OOO; 
Sala bert et al.l l2007). The variations in mode widths and heights 
are also functions of frequency. We have modelled the change in 
profile widths and heights in the manner shown in Fig. [5] Note 
that the size of the variation also depends linearly on the level of 
activity and the profiles shown in Fig.[5]represent the maximum 
changes in width and height respectively, which occur when the 
solar activity is at a maximum. We have varied the frequency of 
the mode in the same manner as before. We have then produced 
instantaneous profiles of a mode for each day in the 8640-d time 
series and we have used these profiles to produce an integrated 
profile which represents that observed in a long time series. 

The total devil-in-the-detail corrections applied to the 8640- 
d BiSON data are plotted in Fig. [6] These corrections include 
both the effect of the distribution of the activity level and the 
cross-talk between the variations in mode frequency, width, and 
height over the observation period. Comparison with Fig. [3]im- 
plies that the distribution of the activity level accounts for ap- 
proximately half of the total devil-in-the-detail correction. The 
frequency corrections shown in Fig. [6] should be added to the 
observed frequencies, v n j, to obtain frequencies that have been 
corrected for the devil-in-the-detail effect. For the 8640-d data 
set considered here the maximum devil-in-the-detail corrections 
are of the order of 5 times larger than the errors associated with 
the raw fitted frequencies. 
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Fig. 6. Devil-in-the-d etail corrections that were used in 
Broomhal fetail (120091) . 



4. Sun-as-a-star correction 

Sun-as-a-star observations are taken from a perspective where 
the plane of the Sun's rotation axis is nearly perpendicular to 
the line-of-sight. Therefore, only modes where I + m is even 
have a non-negligible visibility in Sun-as-a-star observations. 
Furthermore, estimates of the centroid frequencies are domi- 
nated by the \m\ = I components, which are the most prominent 
components observed in the multiplets. It can, therefore, be dif- 
ficult to estimate the true centroid frequency of a mode. 

The difference between the centroid and fitted frequencies 
depends on the level of activity over the period the observations 
were made. When the solar activity is at a minimum the com- 
ponents are observed to be in a near-symmetrical arrangement 
and so the fitted centroid frequency is close to the true centroid 
frequency. However, at moderate to high activity levels this is 
not the case as the mode components are not arranged symmet- 
rically. Take, for example, an / = 2 mode. The \m\ = 2 compo- 
nents will experience a larger shift at high activity levels than the 
m — component. The magnitude of observed asymmetry is re- 
lated to the inhomogeneous distribution of the solar activity over 
the surface and the spherical harmonic associated with each vis- 
ible m component. Therefore, the fitted frequencies differ from 
the true centroids by an amount that is dependent on /. 

lAppourchaux & Chaplin (2007) describe how to make a 
'Sun-as-a-star' correction, which is determined using the so 
called a coefficients that are found from fits for an unre- 
solved Sun. The a coefficients were obtained from MDI data 
dSchou et al.ll l998|). which had th e same activity levels as the 
BiSON data dChaplin et"aIH 2~004a). It is important to make this 
correction as the centroid frequency contains information on 
the spherically symmetric component of the internal structure 
and so is the required input for hydrostatic structure inversions. 
Furthermore, resolved solar observations are able to offer direct 
estimates of the centroid frequencies and so, to compare the fre- 
quencies determined using resolved and Sun-as-a-star data, the 
Sun-as-a-star correction should be employed. 

Figure [7] shows the Sun-as-a-star corrections applied to the 
8640-d BiSON data set. As can be seen, except for 1 = modes 
these corrections are similar in magnitude but opposite in sign 
to the devil-in-the-detail corrections. As the variation with fre- 
quency is different for the two different corrections it is still nec- 
essary to use both when examining the frequencies observed in 
long time series. The correction to the / = modes is zero as 
there is only one component to observe and by definition this is 
the centroid frequency. 
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Fig. 7. Sun-as-a-star corrections that were used in 
Broomhal l et all d2009t) . 



5. Summary 

All three of the corrections described in this paper are depen- 
dent on the solar activity level over the period of the observa- 
tions and so should be used when comparing the frequencies 
determined from data observed at different epochs. The linear 
solar cycle correction allows frequencies to be determined for 
a nominal activity level. The devil-in-the-detail correction ac- 
counts for the distribution of activity levels over the period of 
observations and the fact that mode profile heights, widths, and 
frequencies all vary with solar activity. This correction should be 
used when comparing frequencies obtained from data sets of dif- 
ferent lengths, although the correction is only significant when 
the data sets are long. The Sun-as-a-star correction allows for 
the fact that only components where / + m is even can be de- 
tected in Sun-as-a-star data and should be used when comparing 
the frequencies obtained from resolved and Sun-as-a-star data. 
This is because, at times of high-activity, the m-components of a 
mode are not arranged symmetrically, and so, because not all 
components are visible, Sun-as-a-star observations measure a 
frequency that is slightly different to the centroid frequency. 
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